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ABSTRACT 

The diffuse gamma-rays in the Galactic center region have been studied. We propose 
that there exists a population of millisecond pulsars (MSPs) in the Galactic Cen- 
ter, which will emit GeV gamma-rays through the synchrotron-curvature radiation as 
predicted by outer gap models. These GeV gamma-rays from unresolved millisecond 
pulsars probably contribute to the diffuse gamma-ray spectrum detected by EGRET 
which displays a break at a few GeV. We have used the Monte Carlo method to obtain 
the simulated samples of millisecond pulsars in the Galactic center region covered by 
EGRET (~ 1.5°) according to the different period and magnetic field distributions 
from the observed millisecond pulsars in the Galactic field and globular clusters, and 
superposed their synchrotron-curvature spectra to derive the total GeV flux. Our sim- 
ulated results suggest that there probably exist about 6000 unresolved millisecond 
pulsars in the region of the angular resolution for EGRET, whose emissions could 
contribute significantly to the observed diffuse gamma-rays in the Galactic center. 

Key words: Galaxy: center - gamma-rays: theory - pulsar: general - radiation 
mechanisms: nonthermal 



1 INTRODUCTION 



Gamma-ray emissions liave been detected from tlie Galac- 
tic center (GC) region (Churazov et al. 1994). EGRET on 
board the Compton GRO has identified a central (< 1°) 
~ 30MeV - 10 GeV continuum source (2EG J1746-2852) 
with a luminosity of ~ 10'''^ erg s^^ (Mattox et al. 1996). 
Further analysis of the EGRET data obtained the diffuse 
gamma ray spectrum in the Galactic center. Allowing for a 
total source-excess extent up to 1.5 degree in radius, the 
luminosity (> 100 MeV) attributed to the source excess 
at GC is about 2 x lO^'^erg s~^ (Mayer-Hasselwander et al. 
1998). The photon spectrum can be well represented by a 
broken power law with a break energy at ~ 2 GeV. Be- 
low this energy the diff'erential photon spectrum is F{E) = 
2.2 X 10"^°(£;/1900MeV)"^-^, above the break energy the 
spectrum is F{E) = 2.2 x 10"^°(£;/1900MeV)"^ ^ A re- 
analysis of EGRET data by Hooper and Dingus (2002) has 
shown that the EGRET GeV source is displaced from the 
GC. Recently, Tsuchiya et al. (2004) have detected sub- 
TeV gamma-ray emission from the direction of the Galac- 
tic Center using the CANGAROO-II Imaging Atmospheric 
Cherenkov Telescope. Their data suggest that the GeV 
source 3EG 1746-2851 may be coincident with this TeV 
source. Recent observations of the GC with the air Cerenkov 
telescope HESS (Aharonian et al. 2004) have shown a signif- 



icant source centered on Sgr A* above energies of 165 GeV 
with a spectral index F = 2.21 ± 0.09. 

Some researchers have studied the possible origin of the 
gamma-rays from GC. Mastichiadis & Ozernoy (1994) pro- 
vided that the gamma-rays originate close to the massive 
black hole (Mbh ~ W^Mq), possibly from relativistic par- 
ticles accelerated by a shock in the accreting plasma. In 
the same time, the gamma-rays could come from some ex- 
tended features like radio arcs, where relativistic particles 
are present (Pohl 1997). Markoff et al. (1997) discussed in 
detail the gamma-ray spectrum of GC produced by syn- 
chrotron, inverse Compton scattering, and mesonic decay 
resulting from the interaction of relativistic protons with hy- 
drogen accreting onto a point-like sources (e.g. the massive 
black hole). 

However, the above models cannot produce the hard 
gamma-ray spectrum with a sharp turnover at a few GeV, 
which is observed for the GC source. On the other hand, 
the spectrum is similar with the gamma-ray spectrum emit- 
ted by middle-aged pulsars and millisecond pulsars (Zhang 
& Cheng 2003; Cheng et al. 2004a). Because of the high 
gamma-ray luminosity observed in GC, we may require a 
pulsar population in the inner region. In the following, we 
will first argue that canonical pulsars (including young and 
mature pulsars) may not be a major contributor to the pul- 
sar population in GC. 

Young pulsars are not likely to be a major contribu- 
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tor since few supernova remnants are presently observed in 
the Galactic center field targeted in the deep X-ray surveys 
(Wang ct al. 2002a; Muno ct al. 2003). This viewpoint is 
also supported by pulsar birth rate estimates. Specifically, 
the birth rate of young pulsars in the Milky Way is about 
1/150 yr (Arzoumanian, Chernoff, & Cordes 2002). As the 
mass in the inner 20 pc of the Galactic center is ~ 10* 
(Launhardt, Zylka, & Mczgcr 2002), the birth rate of young 
pulsars in this region is only 10""^ of that in the entire Milky 
Way, or ~ 1/150 000 yr. We note that the rate may be in- 
creased to as high as ~ 1/15000 yr in this region if the star 
formation rate in the nuclear bulge was higher than in the 
Galactic field over last lO'' - 10* yr (see Pfahl et al. 2002). 
Few young pulsars are likely to remain in the Galactic cen- 
ter region since only a fraction (~ 40%) of young pulsars 
in the low velocity component of the pulsar birth velocity 
distribution (Arzoumanian, Chernoff, & Cordes 2002) would 
remain within the 20 pc region of the Galactic center studied 
by Muno et al. (2003) on timescales of ~ 10^ yrs. Mature 
pulsars can remain active as gamma-ray pulsars up to Iff yr, 
and have the same gamma-ray power with millisecond pul- 
sars (Zhang et al. 2004; Cheng et al. 2004a), but according 
to the birth rate of pulsars in GC, the number of gamma-ray 
mature pulsars is not higher than 10. 

On the other hand, there may exist a population of old 
neutron stars with low space velocities which have not es- 
caped the Galactic center (see Belczynski & Taam 2004). 
Such neutron stars could have been members of binary 
systems and been recycled to millisecond periods, having 
formed from low mass X-ray binaries in which the neu- 
tron stars accreted sufficient matter from either white dwarf, 
evolved main sequence stax or giant donor companions (Bel- 
czynski & Taam 2004, in preparation) . The current popula- 
tion of these millisecond pulsars may either be single (having 
evaporated its companion) or have remained in a binary sys- 
tem. Cheng et al. (2004b) provide that wind nebulae of the 
millisecond pulsar population in the Galactic center region 
can contribute to the unidentified X-ray sources in GC by 
the Chandra survey (Wang et al. 2002a; Muno et al. 2003). 
Because millisecond pulsars also remain active as gamma- 
ray pulsars, radiating gamma-rays through the synchrotron- 
curvature process (Cheng & Zhang; 1996; Zhang & Cheng 
2003), it is possible that the observed gamma-ray luminos- 
ity in GC may be produced through an accumulation of 
these millisecond pulsars which would provide the observed 
gamma-ray spectrum. 

In this paper, we will examine in detail if the millisecond 
pulsar population could contribute to the diffuse gamma-ray 
spectrum iu the Galactic center region. In §2, we will present 
our motivation to consider the contribution of millisecond 
pulsars to diffuse gamma-ray in GC based on the outer gap 
model. To find the gamma-ray spectrum of these millisecond 
pulsars, we assume their globular parameters like the period, 
magnetic field are similar to those of the observed millisec- 
ond pulsars. In §3, we derive the period, magnetic field dis- 
tributions of total MSPs, MSPs in the Galactic field and in 
globular clusters from the present pulsar survey database. 
In §4, we sample millisecond pulsars by the Monte Carlo 
method according to the different distributions, superpose 
their spectral profiles to fit the observed diffuse gamma-ray 
spectrum in GC and obtain the number of MSPs which are 
needed in the Galactic center region covered by EGRET. 



Our results are summarized and discussions are also pre- 
sented in §5. 



2 MOTIVATIONS 

Since the mass within the region with the radius ~ 20 
pc (17' X 17') of the Galactic center is ~ 10* M© (Laun- 
hardt, Zylka & Mezger 2002), an estimate for the fraction 
of millisecond pulsars in this region is about 10"'^ of the 
entire Galaxy. Based on the population analysis of Lyne et 
al. (1998), the number of millisecond pulsars in the entire 
Galaxy may exceed 2 x 10®, suggesting that more than 200 
millisecond pulsars exist in the Galactic center region if their 
evolutionary formation channels are similar to the rest of the 
Galaxy. Coupled with the fact that the escape velocity from 
the Galactic center is about 200km s~^ and the average birth 
velocities of observed millisecond pulsars are ~ 130km s~^ 
(Lyne et al. 1998), these pulsars are likely to remain in the 
Galactic center through their entire lifetime. Belczynski & 
Taam (2004) have considered the binary population syn- 
thesis in the Galactic center region, and their results show 
that there exist about 300 low-mass binary systems in the 
Galactic center. Furthermore, about 100 - 200 millisecond 
pulsars could be produced through the recycle scenario and 
lie in the Galactic center region (17' x 17', Taam 2004, pri- 
vate communication). Recently, Pfahl & Loeb (2004) pro- 
pose that ~ 1000 radio pulsars may presently orbit Sgr A* 
with periods of < 100 years, in which 1-10 may be detected 
by current radio telescopes. Therefore, we believe that there 
should exist a millisecond pulsar population in the Galactic 
center, which can contribute to the high energy emissions, 
e.g. X-rays, gamma-rays which are detectable. 

Millisecond pulsars can remain active as gamma-ray 
pulsars through their lifetime according to outer gap mod- 
els which are originally proposed by Cheng, Ho & Ruderman 
(1986). Based on the model, Zhang & Cheng (1997) have de- 
veloped a self-consistent mechanism to describe the high en- 
ergy radiation from spin-powered pulsars. In the model, rel- 
ativistic charged particles from a thick outer magnetospheric 
accelerator (outer gap) radiate through the synchrotron- 
curvature radiation mechanism (Cheng & Zhang 1996) 
rather than the synchrotron and curvature mechanisms in 
general, producing non-thermal photons from the primary 
pairs along the curved magnetic field lines in the outer 
gap. The characteristic energy of high energy photons emit- 
ted from the outer gap is determined by the global pul- 
sar parameters, including the spin period P, the dipolar 
magnetic field B, and the fractional size of the outer gap 
/ ~ ^.hP'^^/'^^B'^''' (Zhang & Cheng 1997), which is the 
ratio between the mean vertical separation of the outer gap 
boundaries in the plane of the rotation axis and the magnetic 
axis to the light cylinder radius. Then the characteristic 
synchrotron-curvature emission energy is given by (Zhang 
& Cheng 1997) 

~ 5 X 107^/^S?/^p-^/^(-^)-"/*eV, (1) 
Rl 

where B12 is the dipolar magnetic field in units of 10^^ G, 
Rl = cP/2'n is the light cylinder radius, and r is the distance 
to the neutron star. The gamma^ray spectrum drops expo- 
nentially beyond the energy E.^. This self-consistent model 
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has also been developed to describe gamma-ray emission 
from millisecond pulsars (Zhang & Cheng 2003). 

Zhang & Cheng (1998) have studied the contribution 
to the Galactic diflFuse gamma-rays from the unresolved 
spin-powered pulsars using the outer gap model. Their re- 
sults show that the gamma-ray emission from these pulsars 
could contribute significantly to the observed Galactic dif- 
fuse gamma-ray spectrum above 1 GeV. Therefore, we be- 
lieve that a laxge number of millisecond pulsars which lie in 
the Galactic center region could also contribute significantly 
to the diffuse gamma-ray spectrum from the Galactic cen- 
ter. Furthermore, according to the results of Zhang & Cheng 
(2003), the gamma-ray spectral cut-off at ^ a few GeV is 
consistent with the observed spectral properties of diffuse 
gamma-rays in the Galactic center. 

To study the contribution of millisecond pulsars to the 
diffuse gamma-ray radiation from the Galactic center in de- 
tail, e.g. fitting the spectral properties and total luminos- 
ity, we firstly need to derive the period and surface mag- 
netic field distribution functions of the millisecond pulsars 
respectively. And then we integrate contributions from all 
the millisecond pulsars with different periods and surface 
magnetic fields to derive the predicted difi'use gamma-ray 
spectrum, which can be compared with the observed spec- 
trum in the Galactic center region to calculate how many 
MSPs are needed. This is the aim of the present paper. 



3 DISTRIBUTION FUNCTIONS OF 
MILLISECOND PULSARS 

We obtain the period and surface magnetic field distribu- 
tion functions of millisecond pulsars from the observed pul- 
sar data. Here, millisecond pulsars are defined as the pulsar 
with P < 10 ms, and Bsurf < lO"' G. So we find 86 detected 
millisecond pulsars from the latest ATNF Pulsar Catalog ^ , 
in which 41 millisecond pulsars are in the Galactic field, and 
45 pulsars in globular clusters. Since the millisecond pulsars 
in the Galactic field and globular clusters may have different 
properties, we derive the distribution functions of total mil- 
lisecond pulsars at first, and then the distribution functions 
of the millisecond pulsars in the Galactic field and globular 
clusters separately. We should notice that the period deriva- 
tive measurement of millisecond pulsars in globular clusters 
is quite difficult and uncertain, so many of them have not 
been given the surface magnetic field, and the magnetic field 
varies from 10® G to 10^° G. On the other hand, the mea- 
surement of the millisecond pulsars in the Galactic field is 
relatively accurate and reliable, the derived surface magnetic 
field is lower than 10® G. 

When fitting the normalized distribution profile, we 
take the Gaussian function as the following form 



f{x) =fo + 



A 



:exp[-2(: 



W 



(2) 



To derive the period distribution function, we define x — P 
and 1ms < P < 10 ms, and for surface magnetic field dis- 
tribution function, x = logSgurf. If including the total mil- 
lisecond pulsars, the fitting parameters for the period dis- 



tribution are /o = 0.045,^ = 0.60, W = 1.95, = 3.91; 
for magnetic field distribution, the parameters are fo = 
0.015,^ = 0.17, W = 0.57, = 8.49 (7.8 < logP < 10). 
The fitting profiles of the pulsar data have been shown 
in Figure 1. Just including the millisecond pulsars in the 
Galactic field, the parameters for the period distribution 
are fo = 0.059, A = 0.48, W = 2.34, Xc = 4.28; the param- 
eters for magnetic field distribution are /o — —0.06, A — 
0.24, W = 0.65, Xc = 8.40 (7.8 < logS < 9) (Figure 2). Only 
including the millisecond pulsars in globular clusters, the pa- 
rameters for the period distribution are fo = 0.042, A — 
0.67, W = 1.62, Xc = 3.71; the parameters for magnetic 
field distribution are fo = 0.04, A = 0.10, = 0.37, Xc = 
8.71 (7.8 < logB < 10) (Figure 3). 



4 SPECTRAL MODELLING OF DIFFUSE 
GAMMA-RAYS OF THE GALACTIC 
CENTER 

As discussed in Section 2, there exists the population of 
millisecond pulsars in the Galactic center region. Firstly, we 
assume the number of MSPs, A^msp, in GC within the angu- 
lar resolution size of EGRET ^ 1.5°, each of them with an 
emission solid angle Af2 ~ 1 sr and the 7-ray beam pointing 
in the direction of the Earth. We performed calculations to 
sample the parameters (period and magnetic filed) of these 
MSPs by the Monte Carlo method using the above three dis- 
tributions of the observed MSPs derived in Section 3. Here, 
we have assumed no evolution for these millisecond pulsars, 
and they will remain to lie in the Galactic center because 
of their low average proper motion velocity (see Lyne et al. 
1998; Arzoumanian, Chernoff, & Cordes 2002) 

Zhang & Cheng (2003) have proposed a model to de- 
scribe the X-ray and 7-ray emission from MSPs with outer 
gaps. We first calculate the the fractional size fm of the 
outer gap in our simulated MSPs, if fm < 1, the outer gap 
can exist and then the MSP can emit high energy 7-rays, 
fm can be estimated by (Zhang & Cheng 2003) 



; 7.0 X 10"^P 



,26/21 / 



B 



lO^G 



) ' 0; 



2/7 



(3) 



where 5r is the distance where the local magnetic field equal 
to the dipole field, in the following calculations, we assume 
= 5r/W^ ~ 1. The 7-rays are produced in the outer 
gap by synchrotron-curvature radiation, the 7-ray differen- 
tial flux observed on the Earth of the ith MSP can be cal- 
culated by 



AD.CP dE-ydt 



(4) 



where d is the distance of the Galactic center to us, taken as 
8.0 kpc, the solid angle of 7-ray beam Af2 ~ 1 sr, the spec- 
trum d^N/dE-ydt can be calculated according to the equa- 
tion (57) of Zhang & Cheng (1997). The total flux which 
contributes to the 7-ray emission from the GC region can 
be obtained by superposing all MSPs with outer gaps 



(5) 



^ http://www.atnf.csiro.au/research/pulsax/psrcat/ 



where n is the number of simulated MSPs with outer gaps. 
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During our calculations, we let the number of millisec- 
ond pulsars iVivisp as a free parameter to fit the observed 
data points using three different distributions of the pe- 
riod and magnetic field separately. We find that about 6000 
MSPs could significantly contribute to the observed GeV 
flux in the Galactic center region. The calculated proflles 
of superposed spectra of the all millisecond pulsars with 
outer gaps in the Galactic center region according to three 
different distributions of the period and magnetic field are 
shown in Figure 4. The solid line corresponds to the distri- 
butions derived from the total detected millisecond pulsars; 
the dashed line just includes the millisecond pulsars in the 
Galactic field; and the dotted line just includes the millisec- 
ond pulsars in globular clusters. Our predicted spectra are 
consistent with the observed results which have been anar 
lyzed by Mayer-Hasselwander et al. (1998) and Hartman et 
al. (1999). 

In Figure 4, one can find that the predicted spectra 
calculated due to the distributions derived from the total 
observed millisecond pulsars and those just in globular clus- 
ters fit the observed data better. The reduced values of 
three curves fitting to the eight data points are 1.51 (solid 
line), 3.95 (dashed) and 1.62 (dotted) respectively. Our re- 
sults probably imply that the unresolved millisecond pulsars 
in GC will follow the period and magnetic field distributions 
similar to the forms observed in globular clusters. However, 
the number of discovered millisecond pulsars is very lim- 
ited at present, so the statistics may be quite uncertain. 
In addition, the predicted flux is dependent on some pul- 
sar global parameters, so it is too early to make the con- 
clusion. Meanwhile, it should be pointed that the different 
profiles of predicted spectra could also be induced by the 
different number of simulated millisecond pulsars with outer 
gaps (satisfying the criterion fm < 1): about 4000 simulated 
MSPs have outer gaps according to the distributions just 
including the MSPs in the Galactic fleld; while, about 5000 
MSPs can have outer gaps for the derived distributions of 
total observed MSPs and those in globulax clusters. For the 
distributions of MSPs in the fleld, even if we increase the 
number to A'msp ~ 10000, the predicted spectrum fits the 
data not better than the solid and dotted lines in Figure 4. 
After checking the simulated data, we find that according 
to the distributions in the field, the number fraction of sim- 
ulated MSPs which satisfies the criterion < 1 is small, 
and the average 7-ray luminosity of MSPs is relatively low. 

Furthermore, the fraction size of the outer gap /m will 
determine the total flux of millisecond pulsars, i.e. Fj oc f^. 
Zhang et al. (2004) considered the effect of the magnetic 
inclination angle in the calculation of fm, and found that 
fm for the large magnetic inclination angle can increase to 
2-3 times of the original value. Millisecond pulsars generally 
have a large magnetic inclination angle (Ruderman 1991), 
then they can emit the higher gamma-ray flux, so the re- 
quired number of millisecond pulsars in the region with the 
radius of ~ 1.5° to produce the diffuse gamma-rays in GC 
can decrease significantly. Even considering part of millisec- 
ond pulsars have no outer gaps when fm is larger, the simu- 
lated number of millisecond pulsars which can fit the spectra 
best is iVivisp ~ 1000 for the distributions of MSPs in globu- 
lar clusters and total observed MSPs, A^msp ~ 2000 for the 
distributions of MSPs in the field. Probably, these decreas- 



ing numbers of MSPs in GG axe more reasonable both for 
the theoretical prediction and observations. 



5 DISCUSSIONS AND CONCLUSION 

In the present paper, we have studied the diffuse gamma- 
rays in the Galactic center region detected by EGRET. We 
propose that there exists the population of millisecond pul- 
sars in the Galactic center, which will emit gamma-rays 
through synchrotron-curvature radiation predicted in the 
outer gap models. The gamma-ray spectrum of millisecond 
pulsars shows a break at a few GeV, which is similar to the 
spectrum of diffuse gammarrays in GG. Therefore, these un- 
resolved millisecond pulsars could contribute to a significant 
fraction of diffuse gamma-rays in GC. 

According to the period and magnetic field distributions 
of the observed millisecond pulsars, we sampled the global 
parameters (period and magnetic field) of the millisecond 
pulsar population in the Galactic center region by the Monte 
Carlo method. Since the millisecond pulsars in the Galac- 
tic field and globular clusters may have different properties, 
we find three classes of distributions which include total ob- 
served MSPs, just include the MSPs in the field and those 
in globular clusters, respectively. Then we used three possi- 
ble MSP samples to calculate their gamma-ray differential 
spectra, and superpose the profiles to fit the observed data. 
The modelled results suggest about 6000 MSPs are needed 
to match the observed gamma-ray spectrum. We also find 
the superposed spectra of MSPs could fit the observed spec- 
trum well except for the sample derived from the distribu- 
tions of the MSPs in the field (see Figure 4), which probably 
suggests the unresolved millisecond pulsars in the Galactic 
center follow the distributions of the period and magnetic 
field similar to the forms in globular clusters. However, be- 
cause the number of millisecond pulsars in globulax clusters 
with the period and period derivative measurements is lim- 
ited, and the number of millisecond pulsars in GC is also 
unknown, we cannot conclude that the hypothetical millisec- 
ond population in GC could resemble millisecond pulsars in 
globular clusters at present. Furthermore, if the effect of the 
magnetic inclination angle is considered in the calculation 
of the fraction size of the outer gap fm, at least about 1000 
millisecond pulsars for the different distribution are still re- 
quired in the region of the radius ~ 1.5° to fit the diffuse 
gamma-ray spectrum in the Galactic center. 

The multiwavelength observations have shown the com- 
plex structure in the GC region (e.g. Purcell et al. 1997; 
Mayer-Hasselwander et al. 1998; Wang et al. 2002a; Maeda 
et al. 2003) , so different scenarios for the origin of the diffuse 
gamma-rays have been considered as mentioned in §1 (also 
see Mayer-Hasselwander et al. 1998 and references therein). 
Recently, Fatuzzo & Melia (2003) have attributed the GeV 
emission to tt" decay resulting from high energy protons in- 
teracting with the ambient matter in Sgr A East. As pointed 
out in Mayer-Hasselwander et al. (1998), those models can- 
not easily produce the very hard spectrum with a sharp 
turnover. In this paper, we suggested that the spectrum 
turnover may be contributed by the millisecond pulsar pop- 
ulation in GC through curvature-synchrotron radiation in 
the magnetosphere. The predicted spectrum of millisecond 
pulsars show a cutoff above ~ 3 GeV (see Figure 4), so 



© 0000 RAS, MNRAS 000, 000-000 



Contribution to diffuse gamma-rays in the Galactic center region from unresolved millisecond pulsars 5 



it cannot contribute to the sub-TeV flux recently detected 
by CANGAROO-II (Tsuchiya et al. 2004) and HESS (Aha- 
ronian ct al. 2004), LtcV ~ lO^'^erg s"^ Those ToV pho- 
tons in GC are possibly induced by pion decay produced 
through pp interactions (e.g. Tsuchiya et al. 2004 and ref- 
erences therein), and compact wind nebulae of millisecond 
pulsars through inverse Compton scattering (e.g. Aharonian, 
Atoyan & Kifunc 1997; Wang ct al. 2004). According to the 
estimation of Wang et al. (2004), the TeV luminosity of the 
MSP wind nebula through inverse Compton scattering is 
about lO'^^erg s^^, with < Lsd >~ 10^'*erg s~^, the aver- 
age medium density n --^ lO^cm^'^ and magnetic field in 
the Galactic center region B ~ 50/iG (Uchida & Giisten 
1995), then the total TeV luminosity contributed by wind 
nebulae of MSPs is 6 x lO^^'erg s~^. Compared with the 
present observations, wc think compact wind nebulae of mil- 
lisecond pulsars through inverse Compton scattering could 
contribute to the TeV flux in the Galactic center, and the 
photon index V ~ 2.2 is also well within the predicted range 
by wind nebula models F ~ 2 — 2.5 (e.g. Wang et al. 2004). 
Hence, the origin of high energy gamma-rays (GcV - TeV 
energy band) is still a mystery, requiring the further obser- 
vational constraints. 

Michelson et al. (1994) have derived an upper limit to 
> 100 MeV luminosity of the globular cluster 47 Tuc, where 
over 20 millisecond pulsars have been idcritifled and an esti- 
mated total population could be larger than 200 (Camilo et 
al. 2000). Using the EGRET observed upper limit of 47 Tuc 
(~ 1.2 X 10^'"' erg s^"*^ together with the estimate MSP pop- 
ulation, the estimated gamma-ray luminosity for individual 
MSP is roughly 6 x lO^'^ erg s~^, which is lower than our 
estimate in this paper. However, Cheng and Taam (2003) 
have studied the X-ray properties of the MSPs in 47 Tuc, 
in which the average X-ray luminosity per MSP is about 
factor 10 lower than those typical MSPs in the held and 
the spectrum is dominated by thermal spectrum with an 
unusually high polar cap temperature. They conclude that 
all these unusually properties can be explained by the fact 
that strong, small-scale, multipole magnetic fields exist on 
the surface of MSPs in the 47 Tuc. They suggest that the 
typical age of MSPs in globular cluster is much older than 
those in the field, therefore the surface magnetic field buried 
under the crust during the accretion phase may difi^use back 
to the surface in billion year time scales. They also suggest 
that these complicated surface magnetic field structure can 
quench the outer gap because gamma-rays emitted from the 
polar gap can become pairs in those magnetic field lines con- 
nected with outer gap (Ruderman and Cheng 1988). There- 
fore they suggest that MSPs in 47 Tuc are weak gamma-ray 
emitters. 

Here, we would like to emphasize two important points 
again. First, the evolution of magnetic field of millisecond 
pulsars in 47 Tuc could be very special. In section 5 of Cheng 
and Taam (2003), who have argued that the density of star 
in 47 Tuc is so high that there could be more than one ex- 
change collision or tidal capture to form a close interacting 
binary system in 47 Tuc. Consequently, the spin-up (or spin- 
down) evolution during accretion (or post accretion) phase 
of a millisecond pulsar in 47 Tuc can be very much different 
from the evolution of other millisecond pulsars. They argue 
that the results in complicated surface magnetic field struc- 
ture quench the outer magnetospheric gaps as suggested by 



Ruderman and Cheng (1988). Therefore, the low gamma-ray 
flux from 47 Tuc may not imply all millisecond pulsars in 
globular cluster emitting weak gamma-ray flux. Secondly, if 
we allow number of millisecond pulsars in the held increase 
by a factor of 50%, then the population of millisecond pul- 
sars in the field can give an equally good fit as those two 
previous cases. Therefore we cannot conclude that the hy- 
pothetical millisecond population in GC should resemble ci- 
ther millisecond pulsars in the held or millisecond pulsars in 
globular cluster from the current data. 

Wc would like to remark that the spectral break of the 
gamma-ray spectrum depends on three pulsar parameters, 
i.e. period, magnetic field and the inclination angle {a) of 
pulsars. The flrst two parameters can be determined in very 
good accuracy, however, the last parameter is difficult to 
be determined accurately. Zhang and Cheng (2003) have 
calculated some model dependent gamma-ray spectra for 
MSPs. They have chosen various inclination angles for differ- 
ent MSPs. For PSR J0437-4715 and PSR J2124-3358, they 
choose larger inclination angle {a ~ 40°), the model spectral 
breaks occur at about 3 GeV (cf. Fig. 1 of Zhang & Cheng 
2003). On the other hand, for PSR J0218-f4232 and PSR 
B1821-24 they choose (a ~ 55°), the spectral breaks occur 
at about 10 GeV(cf. Fig. 2 of Zhang & Cheng 2003). Cheng 
et al. (2004a) have also studied that effect of inclination 
angle on the gamma-ray spectrum, indeed, the spectrum is 
quite sensitive to this parameter (cf. Fig. 1 of Cheng et al. 
2004a). In their study they find that pulsars in the galactic 
plane are younger and have larger inclination angle whereas 
pulsars in medium and high galactic latitude are older and 
have smaller inclination angles. They have used a Monte 
Carlo simulation method to show that if the distribution 
of pulsar inclination angle satisfies a random distribution 
(Biggs 1990), then they find that the gamma-raj' spectral 
break of galactic pulsars arc higher than those in high lati- 
tude by a factor of 3. Although there are very few confirmed 
gammarray millisecond pulsars, i.e. only three possible can- 
didates in EGRET catalog with 3(t level (Fierro 1995) and 
one good candidate PSR J0218+4232 (Kuiper et al. 2000), 
we believe that high energy satellites, like Integral, GLAST, 
can very soon provide more information for us about the 
spectral behavior of MSPs. 

Finally, we are aware that the millisecond pulsar popu- 
lation in the Galactic center region is still an assumption at 
present because there arc no MSPs discovered in the inner 
region. We think the assumption would be reasonable if com- 
pared with the observed MSPs in the Galaxy, specially many 
MSPs discovered in the systems of globular clusters, e.g. 47 
Tuc (Grindlay ct al. 2002), and an estimated total rmmber 
over 200 MSPs (Camilo et al. 2000). In theory, the binary 
population synthesis in the Galactic center region also sug- 
gested that there are hundreds of MSPs in GC (Belczynski & 
Taam 2004; Taam 2004, private communication). However, 
because the electron density in the direction of GC is very 
high (Cordes & Lazio 2002), it is difficult to detect millisec- 
ond pulsars by the present radio telescopes.^ X-ray studies 



^ The Parkes telescope has performed a short-time survey cov- 
ering the region of GC, no pulsars have been reported in 
the initial results (Wang, N. & Manchester, D. 2004, private 
communication) . 
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of the sources in GC would probably be a feasible method to 
find millisecond pulsars by Chandra and XMM-Newton. Re- 
cent deep X-ray surveys of the Galactic center have found 
a large number of unidentified faint X-ray sources (Wang 
et al. 2002a; Muno et al 2003). Cheng et al. (2004b) have 
suggested that synchrotron X-ray nebulae around millisec- 
ond pulsars could contribute to a fraction of these sources, 
and the sources with tailed features would be the good can- 
didates. Some bright X-ray tails probably coincident with 
radio filaments have been suggested to be the pulsar wind 
nebulae after the image and spectral analyses (Wang, Lu & 
Lang 2002b; Lu, Wang & Lang 2003; Sakano et al. 2003). In 
a word, though still in dispute, the millisecond population 
in GC is suggestive, which could contribute to faint X-ray 
sources and diffuse gamma-rays in GC as studied in this 
paper. 

We are grateful to the referee for the comments and 
suggestions, and to N. Wang, R. Manchester, R. Taam, and 
T.C. Weekes for the useful discussions. This work is sup- 
ported by a RGC grant of the Hong Kong Government. 
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Figure 1. The distributions of the period (up) and surface mag- 
netic field (bottom) of total detected millisecond pulsars. The 
dashes lines are the fitting curves using the Gaussian function. 
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Figure 2. The distributions of the period (up) and surface mag- 
netic field (bottom) of the detected millisecond pulsars in the 
Galactic field. The dashes lines are the fitting curves using the 
Gaussian function. 



Figure 3. The distributions of the period (up) and surface mag- 
netic field (bottom) of the detected millisecond pulsars in globular 
clusters. The dashes lines are the fitting curves using the Gaussian 
function. 



© 0000 RAS, MNRAS 000, 000-000 



8 Wang, W., Jiang, Z.J. and Cheng, K.S. 




1E-13 I ' ' ' 

10 100 1000 10000 

Energy (MeV) 



Figure 4. Simulated spectra of millisecond pulsars (assumed 
^MSP = 6000) in our model compared with observed gamma- 
ray energy distribution of the Galactic center region by EGRET. 
The data points are taken from the analyses of Hartman et al. 
(1999) and Maycr-Hasselwandcr et al. (1998). The lines represent 
the integrated fluxes from millisecond pulsars according to three 
different distributions (period and magnetic field) of millisecond 
pulsars assumed in the Galactic center, respectively. The solid line 
corresponds to the distributions derived from the total detected 
millisecond pulsars; the dashed line just includes the millisecond 
pulsars in the Galactic field; and the dotted line just includes the 
millisecond pulsars in globular clusters. 
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